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We discuss a model for the study of quark-hadron duality in inclusive electron 
scattering based on solving the Dirac equation numerically for a scalar confining 
linear potential and a vector color Coulomb potential. We qualitatively reproduce the 
features of quark-hadron duality for all potentials considered, and discuss similarities 
and differences to previous models that simplified the situation by treating either 
the quarks or all particles as scalars. We discuss the scaling results for PWIA and 
FSI, and the approach to scaling using the analog of the Callan-Gross relation for 
y-scaling. 
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I. INTRODUCTION 



Quark-hadron duality was first discovered experimentally in inclusive inelastic electron 
scattering by Bloom and Gilman [lj], more than 30 years ago. In the past fewyears, quark- 
hadron duality has generated a lot of interest both on the experimental 0, y| and on the 
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^Duality is a major point 
It is also the basis for 



theoretical side 

in the planned 12 GeV upgrade of CEBAF at Jefferson Lab 
using QCD sum rule s 1 1911 . and plays an important role in the study of semileptonic decays 
of heavy mesons I20I . bJ]. 

Quark-hadron duality implies that in certain kinematic regions, the appropriate average 
of hadronic observables is described by a perturbative quantum chromodynamics (pQCD) 
result. This is of great practical interest, as we are actually able to carry out a perturba- 
tive QCD calculation, in contrast to a full QCD or full hadronic calculation. Surprisingly, 
duality was experimentally shown to hold in inclusive inelastic electron scattering down to 



momentum transfers of Q 2 ~ 0.5 GeV 2 [2J. Duality also holds in the semi-leptonic decay of 
heavy quarks (2^], and in the annihilation e + e~ —>■ hadrons. The exact manner of averaging 
depends on the process. 

Duality is not only a very interesting phenomenon by itself, but it also has extremely 
important applications. As duality connects the resonance region, i.e. the region where 
the final state invariant mass W < 2 GeV, and the deep inelastic region, one may infer 
information on one from the other. The earliest example discussed was the extraction of the 
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elastic nucleon form factor from the deep inelastic scaling curve 
contributions were inferred from the resonance data. This connection afforded by duality 
opens up the large xsj regime experimentally, as higher Q 2 measurements are difficult to 
obtain - note that the (e, e') cross section contains the Mott cross section as a factor, and 
the Mott cross section is proportional to 1/Q A - The large Xsj region is much easier to access 
in the resonance region, as the necessary Q 2 values there are much smaller. 

One of the most exciting and promising applications of duality will be the measurement 
of the neutron polarization asymmetry A™ at large xbj- There are many different theo- 
retical predictions for this quantity, ranging from (unbroken SU(6)) to 1 (pQCD) 25]. 
Experimental information on A™ at large Xsj would greatly enhance our understanding of 
the valence quark spin distribution functions. There are very recent new data from Jeffer- 
son Lab, going up to Xbj = 0.6 j^J, with small error bars, but the deep valence region of 
Xsj — * 1 remains inaccessible in deep inelastic scattering. If duality is well understood, one 
may take data in the resonance region, apply a proper averaging procedure, and thus obtain 
results for A\{xBj — > 1). 

New theoretical approaches to a better understanding of duality have been based on mod- 
eling: one branch uses the non-relativistic constituent quark model, with some relativistic 



corrections, to describe duality 
a relativistic one-body equation 



Jj. a nd another branch starts the modeling with 

act w th 

the phenomenology. It was started by the pioneering work of Close and Isgur j^J, where 



the authors investigated how a summation over the appropriate sets of nucleon resonances 
leads to parton model results for the structure function ratios in the SU(6) symmetric quark 
model. This work was recently expanded I6j to include the effects of SU(6) spin-flavor sym- 
metry breaking. In the authors considered the first five low-lying resonances. Our 
results belong to the latter branch. The goal of these modeling efforts is obvious: to gain an 
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understanding of quark-hadron duality and the conditions under which it holds, by captur- 
ing just the essential physical conditions of this rather complex phenomenon. We imposed 
these basic requirements for a model: we require a relativistic description of confined valence 
quarks, and we treat the hadrons in the infinitely narrow resonance approximation. 

This paper is the third in a series on modeling quark-hadron duality in inclusive electron 
scattering - the reaction in which quark-hadron duality was first observed by Bloom and 
Gilman. Previously, we simplified the situation by first assuming that all particles involved 
are scalars Q], and then relaxed these constraints for beam electrons and exchange photons, 
and only assumed scalar quarks lfil ]. While these simplifications are physically significant, 
they allowed us to calculate all interesting quantities analytically or semi-analytically. We 
found that the features of quark-hadron duality were reproduced qualitatively in both mod- 
els. Now, we have taken one more step towards a realistic description of the problem: 
previously, we simplified the problem by discussing scalar quarks, but now, we use proper 
spin-1/2 quarks in our model. This lays the foundation for investigating duality in polariza- 
tion observables. 

Also, for the first time, in this paper we present calculations for three different confining 
potentials. Previously, we used a linear potential in the Klein-Gordon equation, which leads 
to a "relativistic oscillator" . Now, we present a scalar linear confining potential and combine 
it with a vector potential: either with a static color Coulomb potential or with a running 
color Coulomb potential. 

The paper is organized as follows: first, we introduce the model and then give model 
results in PWIA. The next section discusses y-scaling: the connection between y-scaling 
and Bjorken scaling, y-scaling results from our two previous models, y-scaling in PWIA and 
FSI, and the sensitivity to different potentials. Then, we discuss our results for sum rules 
in PWIA and FSI. In the next section, we derive the analog of the Callan-Gross relation for 
y-scaling, and investigate the onset of scaling through this relation. Then, we summarize 
our results and give a brief outlook. 



4 



II. THE MODEL 

Our model consists of a constituent quark bound to an infinitely heavy di-quark and is 
represented by the Dirac hamiltonian 

H = a-p + (3(m + V s (r)) + V v (r), (f) 

where the scalar potential is a linear confining potential given by 

V s {r) = br, 6 = 0.18GeV 2 . (2) 

We have used the constituent quark mass in this paper, as our main interest is the study 
of quark-hadron duality, which sets in at rather low Q 2 , experimentally Q 2 ~ 0.5 GeV 2 is 
enough. In this kinematic region, the appropriate degree of freedom is the constituent quark, 
which has acquired mass through spontaneous chiral symmetry breaking. We have used a 
value for the quark mass of m = 258.46 MeV - obtained previously in a fit to heavy mesons 
27| . Calculations will be presented where the vector color Coulomb potential is absent, that 
is V c (r) =0, where the vector potential is the simple static Coulomb potential 

Voir) = (3) 

with a s = 0.181 and where the color Coulomb potential is corrected to allow for the running 
coupling constant in a manner similar to that used by Godfrey and Isgur 28]. This potential 
has the form 
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We assume that only the light quark carries a charge, and we choose unit charge for 
the light quark for simplicity. The inclusive cross section is given by the usual Rosenbluth 
equation 

d 2 ° f Q 4 „ , v / Q 2 2 0\ „ , A 

~ ^Mott < —rRdq, v) + \7n + tan o R t(1^ v ) \ > ( 6 ) 



dtldE moit [ g 4 ' \2q 2 2, 

where CMott is the Mott cross section, q is the three-momentum transfer from the electron 
to the target, v is the energy transfer and Q 2 = q 2 — v 2 . The longitudinal and transverse 
response functions for the model are given by 

R L (q, v) = J2\(*f\ e iqr \* )\ 2 6(u + Eo- E f ) (7) 
/ 

and 

2 

R T (q, i/) = X)EK*/l \*o)\ 2 8(v + Eo- E f ) , (8) 

/ i=l 

where E'o is the ground state energy. In terms of these response functions the structure 
functions are 

W 1 (Q 2 ,u) = ^R T (q,u) (9) 



and 



Q 4 r> I , M Q 2 



W 2 (Q 2 , u) = ^jR L (q, v) + ^R T (q, v) . (10) 



The Dirac wave functions and energy eigenvalues are obtained by integrating the Dirac 
equation using the Runge-Kutta-Feldberg technique and solutions are obtained for energies 
up to 12 GeV with the radial quantum number of n = 200 and \k\ < 70. 

In our model, we excite the bound quark from the ground state to higher energy states, 
and do not allow it to decay. Thus, we do not include any particle production in our model, 
and are strictly quantum-mechanical in this sense. We do not have any gluons in our model, 
either, which means that we do not encounter any radiative corrections. Since the response 
functions consist of a sum of delta functions, we choose to smear out the response functions 
by folding with a narrow gaussian for purposes of visualization. The smeared response 
functions are then given by 

R L[T) ^u)= ^'e-^JT— (g/). (11) 

V^e J-oo 

Before presenting numerical results, we would like to remind our readers that, while the 
present model is more realistic than its predecessors, its results should not be compared 
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quantitatively to inclusive electron scattering from a nucleon. Due to the assumption of an 
infinitely heavy antiquark (or diquark) to which the light quark is bound, our calculation 
most resembles inclusive electron scattering from a B-meson, which has never been measured. 
However, the goal of our work is to gain a qualitative understanding of duality, and the 
current simplification is no impediment to this. 



III. THE PLANE WAVE IMPULSE APPROXIMATION 



The analog to pQCD for this model is the plane wave impulse approximation where the 
bound quark is knocked into the continuum by the absorption of the virtual photon. The 
response tensor for this approximation is 

1 f d?p m 



8 J (2tt)3£ p+<j 
xTr 



8{u + E - E p+q ) 



„7 • (p + q) + m 

Y 2 ^ 7 (7 ' Uv{p) + ns ^ ] 



(12) 



where the vector and scalar momentum density distributions n v (p) and n s (p) are defined in 
terms of the ground state wave function 
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with 
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and 
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n v {p) = ( n° v (p), t—Mp) 

1^1 



nS(p) = r^S(p)^(p); 



7T 



n. 



After performing the angular integrals, the response functions can be written as 



lQn 2 q J]y 

<v + E ) 2 + p 2 - q 2 



dppl (is + E )n v (p) + mn s (p) 



+ 



m 



2p 



Kip) 



(13) 
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(16) 
(17) 



(18) 



7 



and 

Rr(q,u) = — — / dppl {u + E )n° v (p) - mn a (p) 
8vr 2 g J\ y \ { 



8ir 2 q J M 

v + E ) 2 — p 2 + q 2 — m 2 (y + E ) 2 — p 2 — q 2 — m 2 
2q 2pq 



<(P) , (19) 



where 



y = ^{y + E ) 2 - m 2 - q . (20) 



IV. Y SCALING 

Duality implies that we see scaling at high enough energy and momentum transfers, and 
that the results in the resonance region, for lower momentum transfers, oscillate around this 
scaling curve |2j • So, the first point that needs to be checked is the onset of scaling in our 
model. For our previous two models, we were able to show analytically that there is scaling, 
and that the scaling curves for PWIA and FSI coincide exactly This result is in contrast to 
0, where a 30% difference was found between the PWIA and FSI scaling curves. In Q|, a 
different one-body equation, a semi-relativistic Hamiltonian of the form H = \Z\pf 2 + \for 
for massless quarks was used, whereas we have used a Klein-Gordon equation in our previous 
modeling, appropriate for spinless "quarks" with mass. Now that we use the Dirac equation, 
we will again have to investigate the important question if the two scaling curves coincide 
or not. The scaling behavior, and a possible violation of "FSI scaling curve = PWIA scaling 
curve" would have implications for our interpretation of deep inelastic scattering data, where 
one usually assumes that FSI can be neglected. The question if and how scaling may arise 
in the presence of FSI has been discussed in the literature before, both for non-relativistic 
QfiQ and —tic B approaches. Here, we a r e interested in dua lity , and need to 
compare the FSI scaling results with the PWIA scaling results. 

Note that there is a certain arbitrariness in defining "scaling" - have we reached scaling 
only once all curves coincide perfectly, or is a minimal shift in the curves when, e.g., doubling 
the momentum, enough? In the following, our usage is that "scaling" is reached when the 
change in the curve is minimal for a very substantial change in the momentum transfer. 
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A. Bjorken scaling and y-scaling 



In our two previous papers on modelling duality [15, 16], we have presented our results 
for the scaling functions at fixed Q 2 as functions of a scaling variable u. In the Bjorken 
limit, Q 2 — > oo, u goes to UBj = ^ x Bj, where xbj = ^rr- is Bjorken's scaling variable, and 
UBj is the appropriately rescaled version of Xsj for the case of a a target with infinite mass 
M. The values of the energy transfer v accessed for various, fixed values of Q 2 are shown in 
Fig. ^ in the left panel. 

We would like to point out that the previously used simplification of modelling the quark 
as a scalar has allowed us to obtain analytic results for the structure functions, giving us 
access to high energy transfers and high Q 2 without any practical, numerical problems. Now, 
however, we have to rely on a solely numerical solution of the Dirac equation, and cannot 
push to arbitrarily high Q 2 and v values anymore. In particular, we find about 28000 energy 
eigenstates, all below 12 GeV energy. While this is an impressive number of states, it is 
clear from Fig. Q that we will not be able to reach the high Q 2 values found necessary for 
scaling in our model with energy transfers of less than 12 GeV. Note that the peak of 
the structure function 1/W2 is localized around u « 2 — 3, so the higher Q 2 values accessible 
at larger u have little practical significance. 

The matrix elements take their simplest form if calculated as functions of q. Therefore, 
we consider y-scaling in this paper, and use the scaling variable as defined in Eq. lffiJj) . 

V = V~(y + E ) 2 - m 2 - q . 



This variable leads to scaling for fixed |q| = q |12(. The kinematics accessed in a y-scaling 
analysis are shown in Figs. ^ and El Fig. ^ bottom panel, shows the values of the energy 
transfer v accessed in the range from y = —2 GeV to y = 2 GeV - this is the range in which 
we have non-negligible contributions to our response functions and structure functions. Fig- EI 
shows that we are now probing a different kinematic region than in w-scaling: we access 
both spacelike and timelike values of the four-momentum transfer Q 2 . At lower, negative 
y values, Q 2 is spacelike and rather large. In this region, we can easily reach large values 
associated with scaling in u or x-type variables. The accessed Q 2 values decrease slowly 
with increasing y, and at a point where y = E + O(^), we reach the photopoint, Q 2 = 0. 
Note that due to the fact that y is independent of q in first approximation, the photopoint 
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y (GeV) 

FIG. 1: The transferred energy v as a function of the scaling variable u for various values of the 
four-momentum transfer Q 2 (top panel), and as a function of the scaling variable y for various 
values of the three-momentum transfer q (bottom panel). 

is reached at the same y value independent of the considered three-momentum transfer q. 
For y > y Q , we probe the timelike region. Note that one also reaches timelike values of Q 2 
for y < —2q — E + O(-); in this region, however, we find zero strength of the responses, 
and it is practically irrelevant for our case. 



10 



50 
40 
30 



% 20 
10 



o 

10 



— i — 1 — i — 1 — 

q = 2 GeV 

q = 4 GeV 

q = 6 GeV 

q = 8 GeV 
q = 10 GeV 




■2<£ 



-0.5 0.5 
y (GeV) 



1.5 



FIG. 2: The transferred four- momentum Q 2 as a function of the scaling variable y for various 
values of the three-momentum transfer q. 

B. Comparison with previous results 

In order to provide a connection between w-scaling and ^/-scaling, we show our results for 
our two previous models with scalar quarks for y-scaling. We start out with the all-scalar 
results discussed in jl^. In we treated beam, exchange particle and quark target as 
scalars. This simplified treatment allows for an analytic solution of the problem. With all 
particles being scalars, only one structure function appears, there is little structure, and 
scaling sets in around Q 2 ~ 15 GeV 2 . Scaling for the PWIA sets in at very low Q 2 of about 
2 GeV 2 . Keep in mind that we scatter from an infinitely heavy target, so that a comparison 
of numerical values for Q 2 for the onset of scaling with actual nucleon data is meaningless. 

Fig. El (top panel) shows the corresponding ^/-scaling plot for the PWIA. One can see 
clearly that for q = 2 GeV, scaling has already set in. Obviously, there are no resonance 
bumps in the PWIA plot, as the final state is a fictitious "free quark". In the bottom 
panel of Fig. El we show the same PWIA scaling function, plotted for fixed four-momentum 
transfer Q 2 versus u. The overall features of the curves, plotted either way, are the same: 
they are smooth and scale quickly. In the w-scaling plot, one sees that scaling has set in at 
Q 2 = 6 GeV 2 , and the changes from the lower value of Q 2 = 2 GeV 2 to Q 2 = 6 GeV 2 are 



11 




FIG. 3: The scaling function S cq of the all-scalar model [13] is plotted versus y for several values 
of the three-momentum transfer q in the top panel, and it is plotted versus u for several values of 
the four-momentum transfer Q 2 in the bottom panel. The results shown have been calculated in 
PWIA. 

tiny. 

Fig. HJ top panel, shows the ^/-scaling plot for the all-scalar model including FSI. Here, 
one can see that scaling does take a while to set in - the lowest value for q that is displayed, 
q = 2 GeV, still shows the typical resonance bump structure. However, the overall shape 
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FIG. 4: The scaling function S cq of the all-scalar model is plotted versus y for several values of 
the three-momentum transfer q in the top panel, and it is plotted versus u for several values of 
the four-momentum transfer Q 2 in the bottom panel. The results shown have been calculated 
including final state interactions (FSI). 

even at low q is already very close to the scaling curve, but slightly shifted towards higher 
y values. The scaling curve is approximated reasonably at q = 6 GeV, and the curves start 
to coincide with each other once q = 10 GeV is reached. 

This is contrasted with the w-scaling plot for the all-scalar model including FSI in the 
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bottom panel of Fig. |3J One can see how the resonance bumps at lower Q 2 give way to 
smooth curves at larger values of the four-momentum transfer. One can see that scaling 
sets in just in the same way, independent of which set of kinematic variables one chooses for 
plotting. For practical reasons - the current matrix elements take their simplest form as a 
function of q - we choose to present our results as functions of q and y. 

Now, we will proceed to show the y scaling results for the "scalar quark model" jl^ . 
too. This will allow us to study the transition from a simple model to a more sophisticated 
version, and to point out differences and common features. In we included the proper 
spins for the beam electrons and the exchange photons, but still treated the quark as a 
scalar. We will refer to this model as the "scalar quark model" in the present paper. This 
model has a conserved current, and a much richer structure, due to the fact that the photon 
can have transverse and longitudinal polarization, which leads to two structure functions. In 
Fig- El top panel, the PWIA results definitely scale more slowly than for the all scalar model. 
Scaling is reached at q = 20 GeV. Not surprisingly, the FSI results, shown in the bottom 
panel of Fig. exhibit an even slower scaling. Even from q = 20 GeV to q = 50 GeV, a 
small difference is visible in the scaling curves. This is the slowest onset of ^/-scaling which 
we have encountered so far, and this mirrors precisely what we expected from our previous 
w-scaling results. Now we are in a good position to progress to the main topic of this paper, 
the results for proper spin-1/2 quarks. Our study of the onset of scaling for PWIA and FSI 
will allow us to draw some conclusions about the onset of u scaling for the FSI case, even 
if we can't calculate for the necessary Q 2 numerically. First, we are going to present results 
for the PWIA for three different potentials. 



C. y-scaling in PWIA 

First, we show our results for PWIA, where we expect scaling to set in at lower momentum 
transfers than for the FSI results. For the FSI results, the scaling curve is generated in our 
model by many overlapping resonances, and one needs to reach kinematics where a sufficient 
number of resonances is accessible to see scaling. 

The scaling variable y is defined by (|20|) and by inverting to find v as a function of y and 
q, the response functions can be written as functions of q and y. For example the PWIA 
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y (GeV) 
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FIG. 5: The scaling function S2cq of the "scalar quark" model is plotted versus y for several 
values of the three-momentum transfer q. The results shown in the top panel have been calculated 
in PWIA, i.e. for the bound-free transition. The results shown in the bottom panel have been 
calculated including final state interactions (FSI), i.e. for the bound-bound transition. 
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response functions (JTSj) and (|19p can be rewritten as 

ry+2q 



I ry+2q ( 

RL{q,y) = 1C 9 / rfpps + g) 2 + m 2 n°(p) +mn s (p) 

y 2 + 2qy+p 2 \ 
+ ^ WJj ( 21 ) 



and 



Rr{q,y) = J dpp^ y / (y + q) 2 + m 2 n° v (p)-mn s (p) 

K(p)] • (22) 



y 2 + 2qy + 2q 2 — p 2 y 2 + 2qy — p 2 
2q 2pq 

In the limit of large q the PWIA response functions become 



UmR L (q,y) = / dpp (n°(p) + ^ <(p) 1 (23) 

i67r- y^i [ p j 



and 



lim R T (q,y) = A / ^ - -KiP)] ■ (24) 

3-»do 8^ J lyl I p J 

These response functions therefore scale in y. Note that the longitudinal response is shifted 
toward positive y while the transverse is shifted toward negative y, since n s v (p) is positive. 
The overall peak heights of the longitudinal and transverse responses in PWIA have roughly 
a 1 : 2 ratio. 

We start by investigating the behavior of the PWIA responses for the linear potential: 
V s = br, V v = 0. In Fig. El bottom panel, we show the longitudinal response Rl for various, 
lower three-momentum transfers q. Here, one sees clearly that both the peak height and 
the width of the peak are reduced significantly by increasing q for q < 10 GeV. The 
peak position also shifts very slightly to lower y values. In Fig. El top panel, we show the 
longitudinal response Rl for various, high three-momentum transfers q. There is a small but 
visible change in going from q = 10 GeV to q = 20 GeV. Scaling is reached for q = 40 GeV, 
and the curves corresponding to even higher momentum transfers coincide. 

In Fig.0 we show the corresponding results for the transverse response function Rt- For 
the low momentum transfers, bottom panel, one sees that increasing q leads to an increase in 
peak height and width, along with a small shift of the peak position towards lower y values. 
The top panel with the high momentum transfers shows that scaling sets in faster for Rt 
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FIG. 6: The longitudinal response function Rl is plotted versus y for several low (bottom panel) 
and high (top panel) values of the three-momentum transfer q. The results shown have been 
calculated in PWIA, using the linear potential. 

than for Rl- Changes are smaller from one q value to the next, and the curves coincide once 
q = 40 GeV. 

It is interesting to consider the onset of scaling when making our modeling more realistic: 
in the transition from the "all scalar" to the "scalar quark" model, we observed a considerable 
slowdown in the onset of scaling, due to the additional structure, and to the more complicated 
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FIG. 7: The transverse response function Rt is plotted versus y for several low (bottom panel) and 
high (top panel) values of the three-momentum transfer q. The results shown have been calculated 
in PWIA, using the linear potential. 

form of the terms. Now, we have added the proper spin to the quarks, but interestingly, 
we see no significant change in the scaling behavior. This is illustrated in Fig. |HJ where 
we show high q results for the scalar quark model and our current model. While the two 
models obviously lead to different scaling curves, the onset of scaling occurs at roughly the 
same momentum transfers in both cases. From this, we have to conclude that the spin of 
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FIG. 8: A comparison of the onset of scaling for the "scalar quark" model and the model with 
quark spin. The transverse response function Rt and the scaling function S2cq of the "scalar quark" 
model are plotted versus y. The results shown have been calculated in PWIA, and have already 
been shown in Figs. 17151 

the quark does not play a significant role in the onset of scaling. 



Now we will show the scaling behavior with FSI included. For reasons of numerical 
feasibility, we restrict ourselves to momentum transfers of q < 10 GeV. 

Figure M shows the longitudinal (top panel) and transverse responses (bottom panel) 
as a function of y for momentum transfers from 2 to 10 GeV with a width of e = 0.02 
GeV for smoothing (recall Eq. (jlljO . As the momentum transfer increases, the average 
response moves to lower y and since the density of states is also increasing the curves 
become increasingly smooth. The longitudinal response approaches the asymptotic result 
from above, while the transverse response approaches it from below. The low q curves 
with the visible resonance bumps oscillate around the smoother curves obtained for higher 
momentum transfer. 

Before comparing the scaling curves obtained in PWIA and with FSI, we need to in- 



D. y-scaling in FSI 
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FIG. 9: The longitudinal response function Rl (top panel) and the transverse response function 
Rt are plotted versus y for several values of the three-momentum transfer q. The results shown 
have been calculated including FSI, using the linear potential. 
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vestigate the influence of our smoothing procedure for the delta-function in energy This 
is an artifact of our model, as we assume that the resonances do not decay. In [16], we 
used a Breit-Wigner type smoothing procedure, and the width F had some influence on 
the numerically obtained scaling curves. However, for the scalar quarks discussed in |lfi| . 
we could find an analytic expression for the scaling curve with FSI. This is not the case 
here, so this matter deserves careful investigation. The PWIA calculation does not suffer 
from this problem, as the knocked-out, "free" quark can have any energy, in contrast to 
the fixed-energy resonance final states included in FSI. Therefore, comparing the FSI and 
PWIA scaling curves is not entirely straightforward. 
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FIG. 10: The transverse response calculated for the linear potential in FSI at q = 10 GeV for 
various values of the Gaussian smoothing parameter e. 

In Fig. [TU1 we show the transverse response Rt at q — 10 GeV, including FSI, calculated 
for various values of the Gaussian smoothing parameter e. We show curves for e = 0.02 GeV, 
the value used for all other plots in this paper, and for e = 0.01, 0.04, 0.08 GeV. While a 
smaller value of e leads to less smoothing and visible resonance bumps even at higher q, the 
overall shape, peak position and width are not changed. For larger values of e, there is no 
visible difference to the original curve with e = 0.02 GeV . Therefore, the dependence on 
the smoothing parameter e is so weak that it will not influence our comparison of the PWIA 
and FSI scaling curves. 
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E. Model dependence 

Calculations for the linear-plus-coulomb potential are shown in Fig. ^2 an d for the 
linear-plus-running potential in Fig. El These figures show features consistent with those 
of the linear potential alone. One can see that the linear potential leads to the highest 
peak value for both and Rl, and has the narrowest width, while the linear-plus- running 
potential has the lowest peak height and the largest width. In all cases, we have not quite 
reached the scaling limit yet, but as the changes from q = 6 GeV to q — 10 GeV are small, 
we are not far away from it, either. The low q results clearly oscillate around the high q 
results, as expected from duality. 

Figure EH shows a comparison of the FSI and PWIA calculations at q = 10 GeV for the 
longitudinal and transverse responses. At this momentum transfer the corresponding FSI 
and PWIA calculations differ by small changes in magnitude and peak position. Since the 
calculations are not yet converged to the scaling limit, it is not clear whether this represents a 
failure to scale to the same limit or is simply a manifestation of differing rates of convergence. 

On a fundamental level, this is important as FSI is usually assumed to be negligible in 
the analysis of deep inelastic scattering data. Recently, however, some authors have pointed 
out that FSI - through gluon exchange between fast, outgoing partons and target spectators 
- can make contributions to the leading twist structure functions at small Xsj 33|- Still, the 
changes we see here are small, and one does not make much of a mistake in neglecting the 
effect of all the FSIs combined. 



V. SUM RULES 

Some of the features of this model can be explored further and more precisely by means of 
energy-weighted sum rules. This will allow us to get a better feeling for the bulk features of 
our results, like peak position and width. Moments that become constant at high momentum 
transfer are another signature of quark-hadron duality (^J, and with the sum rules derived 
in this section we can test if our model behaves in this way. 

By using an integral representation of the energy conserving delta function and replacing 
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FIG. 11: The longitudinal response function Rl (top panel) and the transverse response function 
Rt are plotted versus y for several values of the three-momentum transfer q. The results shown 
have been calculated including FSI, using the linear-plus- Coulomb potential. 
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FIG. 12: The longitudinal response function Rl (top panel) and the transverse response function 
Rt are plotted versus y for several values of the three-momentum transfer q. The results shown 
have been calculated including FSI, using the linear-plus-running potential. 
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eigenenergies with the hamiltonian operator, the response functions can be written as 

/°° tit 
—e ivt < y \ e iH(P,*)t e -i<i- Xe -iH(p, x )t e i q - x ^ o > 
-oo 2n 



-oo 



/°° tit 
— e iut < $ \ e iH(P,*)t e -iH(p+ q , X )t^ o > ^ 
-oo 2vr 



and 



2 f°° tit 

R T (q,u) = ^e lut <^o\e tH{Mt a t e~ iq - x e- iH ^' x)t e iq - x a t \^ > 

i=i 27r 

= J2 ^T eWt < *o|e ifr(ft!B)t a i e-" r (* + * !B )*a i |*o > , (26) 

i=l 27F 

where the momentum shift operator has been recognized and used in the last step for both 
response functions. 

The energy weighted sum rules are defined as integrals of the response functions over 
— oo < v < oo weighted with powers of v. This integral can be simplified by writing 

-i—j e lut (27) 

and then integrating by parts to give 

/OO 
dvv n R L (q, v) 
-OO 

(£) \ n 



(28) 



and 



^Tn(q) = [ dvis n R T (q,is) 

J — OO 

= £ (i^X < * \e iH ®rt t a i e- iH(p + q ' ,Byt ai\*o > 
i=i ^ ' 



(29) 



The momentum-shifted hamiltonian for our model is 

H(p + q, x) = a ■ q + H(p, x) . (30) 

As a result, the sum rules can be reduced to nested commutators involving H, ot - q and aij. 
The first three sum rules for the longitudinal response are 

<W<?) = < V'olV'o >= 1 
SLi(q) = < ipo\a ■ q\ip >= 

S L 2(q) = q 2 (31) 
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The first of these is the Coulomb sum rule and simply reflects the conservation of charge. 
We have assumed that the charge of the bound quark is 1 for simplicity. The second 
of these indicates that the response is roughly antisymmetric about v = 0. The third 
indicates that the width of the distribution is increasing as q. This is consistent with 
the case where the longitudinal response consists of peaks at positive and negative energy 
transfers of roughly equal magnitude and opposite signs that are moving away from the 
origin in opposite directions by a distance proportional to q. The existence of the negative 
energy peak is a consequence of the presence of the negative energy solutions to the Dirac 
equation which are not physical in a simple one-body model. 

The positive energy contribution can be isolated by defining the positive energy moments 

POO 

S L (T)n(q) = / duu n R L(T) {q,u). (32) 

J 

Since the n = moment is no longer normalized to one, we can define the average energy 
transfer as 

Sl(t)i 
Sl(t)o 

and the rms variation in energy transfer as 



(») m = IT^ (33) 



Figures El and El show the positive energy n = moments of the longitudinal and 
transverse response functions. For Slo, the results for the linear, linear-plus-coulomb, and 
linear-plus-running potentials are virtually identical except at q > 10 GeV where the moment 
starts to fall off due to finite range of excitation energies that have been calculated. The 
corresponding calculations in PWIA vary at small q due to variations in the available phase 
space, but must approach | for q —>■ oo. It is clear from this figure that neither the FSI 
calculations nor the PWIA calculations have reached their saturation values at the highest 
values of q that have been calculated here. Although it is plausible that the full results may 
be approaching the PWIA values for large q, it is not possible to determine that this is the 
case based on these calculations. 

Similar results can be seen for the transverse moment shown in Fig. El although there is a 
greater variation among the full calculations at lower values of q. In this case the asymptotic 
value of the moment for the PWIA calculations is 1. 
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FIG. 14: The positive energy n = moment of the longitudinal response function, calculated with 
FSI and in PWIA for the linear, linear-plus-coulomb, and linear-plus-running potentials. 



The asymptotic values of (v) L r T ) f° r ^ ne PWIA are given by 
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lim (z/)j 

q— >oo 



q - E 



12tt 2 



and 



lim (u) T = q-Eo- 



dpp n s v (p) 



dpp n s v (p) . 



(35) 



(36) 



Note that the longitudinal and transverse responses are offset from one another by terms 
that depend upon n s v (p). Since the average value of the energy transfer increases linearly 
with q, it is convenient to plot this moment as q — (y) L . This quantity is shown in Fig. 
[TBI for the three FSI and three PWIA calculations used in the previous figures. Here the 
three FSI calculations have very similar average positions while the average positions of the 
PWIA calculations show large differences. This is the result of the sensitivity of the position 
of the peak of the PWIA response to the difference in energy of the bound state and the 
lowest plane-wave state. In the FSI calculation, both the ground and excited states see 
the same potential which seems to limit the size of the shift in peak position for the three 
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FIG. 15: The positive energy n = moment of the transverse response function, calculated with 
FSI and in PWIA for the linear, linear-plus-coulomb, and linear-plus-running potentials. 

different potentials. For purposes of comparison, it is possible to eliminate this disparity in 
the PWIA calculations by introducing a constant vector potential into the Dirac equation. 
This has the effect of simply shifting the spectrum. As a result we can use this shift to place 
all of the ground state energies at the same value. This shift has no effect on the response 
functions for the FSI calculations, but will correct for the differences in phase space in the 
PWIA. In fact, we already have applied this shift above in the section on y-scaling. 

Figure El shows the rms widths for the various calculations of the longitudinal response. 
The asymptotic value of this width in the PWIA is given by 



lim Aul 



12tt 2 



dpp 4 n® (p) 



144vr 4 



dpp 3 n s v (p) 



(37) 



In all three cases the width of the FSI calculation is close to that of the corresponding PWIA 
calculation. This suggests that the width of the response is determined largely by the width 
of ground state momentum distributions. 

While we have not yet reached convergence to the scaling limit at the momentum transfers 
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FIG. 16: The difference q — (v)l between the average energy transfer (v) L and the momentum 
transfer for the longitudinal response, calculated with FSI and in PWIA for the linear, linear-plus- 
coulomb, and linear-plus-running potentials. 



for which we have calculated, one can clearly see that the moments do flatten out with higher 
momentum transfers, qualitatively agreeing with the observations in duality experiments. 



VI. THE "Y-SCALING CALLAN-GROSS RELATION" 

From a;-scaling in the deep inelastic region, the Callan-Gross relation 2xF\(x) = F^lx) 
is known to hold. The physical significance of the Callan-Gross relation is that one scatters 
off spin 1/2 objects, which have a dominant contribution from the magnetization current, i.e. 
in the transverse part of the structure functions. In the real world, the Callan-Gross relation 
needs to be corrected for radiative effects, and is then observed to hold in the experimental 
data from the deep inelastic region. In our model, there are no radiative corrections, as 
we do not allow the production of new particles. Our model is, so far, entirely quantum- 
mechanical in this sense. Therefore, we expect the Callan-Gross relation to hold as is in x or 
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FIG. 17: The rms width Avl for the longitudinal response, calculated with FSI and in PWIA for 
the linear, linear-plus-coulomb, and linear-plus-running potentials. 

u scaling. As we focus on y scaling, we now derive the analog of the Callan-Gross relation, 
using the fact that both the longitudinal and the transverse response scale for large q. 

The structure functions Wi,W 2 are related to the longitudinal and transverse responses 
Rl and Rt by 



1 



-Ri 



and 



W 2 = %R L + %R T 
g 4 2q z 



(38) 



(39) 



The Callan-Gross relation can be obtained from these equations by re-expressing q and v in 
terms of x and Q 2 , and increasing Q 2 at fixed x. In order to find an analogous relation for 
y scaling, we re-express Q 2 and v in terms of q and y, and then increase q at fixed y. The 
kinematic factors in this limit become 

Q 2 



V- 



2q 2 



E -y 



(40) 
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and 



The two functions, F% = MW\ and F2 = 1/W2, in this limit are: 

F 1 = MWt = -MRt -> -MRt (42) 
2 2 

F 2 = z/W 2 = v%R L + i/^ifr -> (E - y)R T (43) 
g 4 2q z 

Direct comparison of Eqs. (|42j) and @3) yields F 2 = F^Eq - y). Substituting the 
definitions of F\ and F2, one obtains the analog of the Callan-Gross relation for ^/-scaling: 



2W 1 (E -y) = uW 2 . (44) 

We now proceed to check if our numerical results fulfill the Callan-Gross relation, i.e. 
if we have reached the scaling region already. First, we show results for PWIA, where 
we can access arbitrarily high q without any numerical problems. In order to display the 
approach to scaling, we show the ratio of the left-hand-side to the right-hand-side of Eq. 
PH) , tcg = 2W\{Eq — y)/vWv. In Fig. [TBI top panel, we show rcc for the linear potential 
for values of the three momentum ranging from q = 10 GeV to q = 1000 GeV. As expected, 
scaling has set in at these high momentum transfers: the deviations of rcc from 1 are very 
small, at the level of 5% or less for q > 60 GeV. It is interesting to note that for higher y 
values, 1/W2 dominates, while for lower y values, 2W\(Eq — y) is larger. As in the calculation 
of the response functions, we have applied the energy shift to adjust the difference in phase 
space in PWIA. Note that we have omitted the regions y < —lGeV and y > 0.7GeV from 
our graphs, as they do not contain any strength in the responses. The latter is the area 
where (Eq — y) undergoes a sign change, leading to a spike in the ratio rcc that is unrelated 
to the interesting physics - there is hardly any strength in the responses for y > 0.5 GeV. 

For lower q values, see bottom panel of Fig. scaling has clearly not yet set in, and the 
deviations of tcg from unity are large. Comparing these results for rcc with the onset of 
scaling observed in the PWIA response functions Rl and Rt, as discussed in Section HV CI 
we see that rcc gives additional information on the scaling behavior. While the response 
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FIG. 18: The Callan-Gross ratio rcc is plotted versus y for several values of the three-momentum 
transfer q. The results shown have been calculated in PWIA with a linear potential. Note the 
different scales employed in the two panels. 

function plots show that the curves coincide for q = 40 GeV, at the latest, the Callan- 
Gross ratio is more sensitive and still shows slight changes for increasing three-momentum 
transfers. 

Comparing the approach of rcc to 1 for different potentials, Fig. EH shows that linear, 
Coulomb, and running potentials have an almost identical scaling behavior, both at low and 
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FIG. 19: The Callan-Gross ratio rcc is plotted versus y for the linear, Coulomb and running 
potential, for three-momentum transfers q = 2, 10, 20, andWOGeV . The results shown have been 
calculated in PWIA. Note the different scales employed in the two panels. 

high q values. 

In Fig. 1201 we show tqg calculated including FSI for q = 10 GeV, the highest value we 
have attained for FSI calculations. The calculations shown are performed for the different 
potentials considered in this paper, and one can see that, just like for PWIA, the results for 
the Callan-Gross ratio rcc are very similar. 
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FIG. 20: The Callan-Gross ratio rcc is plotted versus y for three-momentum transfer q = 10 GeV. 
The results shown have been calculated in FSI with a linear, Coulomb, and running potential. 

This similarity to the behavior of the PWIA curves leads us to a direct comparison of the 
PWIA and FSI results for the linear potential in Fig. |^] The PWIA and FSI agree very 
nicely. For lower q values, the FSI results obviously lead to oscillatory behavior, but they 
average to the smooth PWIA result. The PWIA and FSI results for the Coulomb potential 
also track each other very closely. For the running potential, the differences between FSI 
and PWIA are a bit more visible, but FSI and PWIA results are still very close. 

While this is no proof, this behavior leads one to the conjecture that scaling for FSI 
should set in at the same q value as for PWIA, and that the rate of convergence in both 
cases is the same. Note that this result does not imply that the individual scaling results 
for PWIA and FSI are the same. Quite on the contrary, as we found that the FSI results at 
q = 10 GeV are not that close to the PWIA scaling curve, and as the PWIA results did not 
change dramatically from q — 10 GeV to higher q values, the conjecture of the same rate of 
convergence suggests that the final scaling results in PWIA and FSI might be different. 
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FIG. 21: The Callan-Gross ratio rcc is plotted versus y for three-momentum transfers q = 2 GeV 
and q = 10 GeV. The results shown have been calculated in PWIA and FSI with a linear potential. 

VII. SUMMARY AND OUTLOOK 



In the present paper, we have expanded our modeling of quark-hadron duality to describe 
a more realistic situation: we now include the spin of the quark. Numerically, this is more 
complicated than the all scalar model and the scalar quark model la], where the 
solutions could be found analytically. However, the inclusion of the quark spin paves the 
way to study spin structure functions, which are one of the most promising areas for the 
practical application of quark-hadron duality. 

We have tracked the changes introduced by making our models more realistic. Most 
notably, adding spin to the quark does not seem to influence the onset of scaling: our 
present model shows about the same scaling behavior as the scalar quark model. This is in 
contrast to the major changes in scaling behavior observed when going from the all scalar 
model to the scalar quark model. 

As before, we have seen that the features of duality observed in experiments, namely 
scaling, oscillation of the resonances around the scaling curve, and flat moments at higher 
momentum transfer, are reproduced qualitatively. We have used a constituent quark mass 
of m = 258.46 MeV for our calculations. This number was taken from a fit to heavy mesons 
. However, nothing hinges on using that particular value: we changed our quark mass to 
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m = 10 MeV, in order to have a value reminiscent of a current quark mass, and repeated 
our calculations. It turns out that, while scaling does set in a little faster, there are no 
qualitative changes in the results. 

In spite of the considerable numerical effort, we have not yet reached full convergence in 
FSI. Therefore, the important question if the FSI scaling curve coincides with the PWIA 
scaling curve could not be answered with certainty yet. At the highest q value we reached, 
q = 10 GeV, we still see small differences between FSI and PWIA curves. Further investi- 
gation of the high q region with different methods will be the subject of a future paper. 

Our findings are different from the results in [8j. There, the author used massless quarks 
and a specific assumption about the potential, namely an identical functional form for the 
scalar potential V s and the vector potential V v . The condition V s = V v leads to a simpli- 
fication in the numerics, as upper and lower components decouple. Results reported in ^] 
indicate a much larger difference between FSI and PWIA results at q = 10 GeV. This may 
be due to the potential chosen there. 

For the first time, we have discussed several different potentials. The combination of a 
scalar linear potential and a vector color Coulomb potential (either with a running coupling 
constant or without) can be considered a fairly realistic approximation to nature. We were 
also able to demonstrate that the observed qualitative features of duality- scaling, low q 
duality, convergence of the moments - persist no matter which potential is used. This hints 
at quark-hadron duality as a fairly general property of inclusive electron scattering. While 
the shape of the response functions is clearly influenced by the ground state momentum 
distribution, and therefore by the potential, the rate of convergence to scaling is unaffected 
by the choice of potential, as demonstrated by studying the validity of the y-scaling Callan- 
Gross relation. 

Duality has been thoroughly explored experimentally in the unpolarized nucleon structure 
functions, and new data on polarization observables are available, too, both from Jefferson 
Lab j^J and Hermes at DESY j^. We will apply our model to the calculation of spin 
structure functions next. 
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